Incident-Energy Dependence of the 
Effective Temperature in Heavy-Ion Collisions 
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I. INTRODUCTION 

Experimental data on transversal momentum distri- 
butions of kaons produced in central Pb+Pb or Au+Au 
collisions show an anomalous dependence on the collision 
energy 0, 0, Q . The effective temperature T* , or slope 
parameter of the transversal- momentum spectra [^J, in- 
creases with energy in AGS and RHIC energy domains. 
However, in the SPS energies the effective temperature 
keeps approximately constant. Recently, it was arg ued@ 
that this behaviour might be caused by a modification of 
the equation of state in the transition region between 
confined and deconfined matter, as suggested by Van 
Hove ||| a long time ago, and could be considered as 
new signal of deconfinement in the SPS energy domain. 

Our main object in the present work is to study the 
behaviour of the effective temperature for K + in several 
energy domains. For thispurpose, we apply the recently 
developed SPheRIO 0, code for hydrodynamics in 
3+1 dimensions, using both Landau- type compact ini- 
tial conditions and spatially more spread ones. For the 
latter we used the average over many events generated by 
NeXus ilia. We show that initial conditions given in 
small volume, like Landau-type ones, are unable to repro- 
duce the effective temperature together with other data 
(multiplicities and rapidity distributions). These quanti- 
ties can be reproduced altogether only when using a large 
initial volume with an appropriate velocity distribution. 

Besides, it seems that the increase of T* in the RHIC 
energy domains is caused mainly by the larger expan- 
sion time in the hadronic phase for higher incidente en- 
ergy, which implies lower freeze-out temperature as y/s 
increases |ll],|I^. It seems that, within our analyses with 
NeXus initial conditions, the RHIC incident energies are 
not enough to produce noticeable increase in the trans- 
verse acceleration during the quark-gluon plasma (QGP) 
phase. 



II. HYDRODYNAMICS IN 3+1 DIMENSIONS 



The hydrodynamical code in 3+1 dimensions used here 
is based on the tec hniq ue called Smoothed Particle Hy- 
drodynamics (SPH) UlEl. This is a method which uses 
Lagrangian coordinates. The fluid is represented by small 
volumes called SPH-particles and the equations for the 
fluid evolution become a system of ordinary differencial 
equations for the SPH-particles, in this representation. 

The numerical code SPheRIO [l^ is a suitable im- 
plementation of this method for the relativistic nuclear 
collisions. An advantage of this technique is the conve- 
nience in the study of problems where the geometry is 
highly irregular, as is the case of non-central relativistic 
collisions. A SPH-particle has attached to it conserved 
quantities; in the present version of SPheRIO, the en- 
tropy and the baryonic number are the quantities which 
are kept constant during the fluid evolution. In the SPH 
representation, the entropy density and the baryonic den- 
sity are parametrized as: 
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where i/f (is?) is the entropy (baryonic number) of the 
i-th SPH-particle, 7 = 1/V(1 - 1% - 1% - t 2 v 2 ), and 
W is the interpolating kernel with width h. Here we 
use the hyperbolic coordinates r = y/(t 2 — z 2 ), x, y and 
77 = 1/2 • ln[(i + z)/(t — z)\ which are convenient for a 
system in rapid longitudinal expansion. The equations 
for the SPH-particles are given by: 
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where oj = (v x ,v y ,T 2 v v ). The pressure, P, and the en- gluons. The thermodynamical quantities are given by 
ergy density, e, are related to the entropy density s and 
the baryonic density rib via equation of state. 
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A. Decoupling criterion 



The distribution of final particles is obtained by using 



Cooper-Frye's prescription 
d 3 N 
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where f(p-u) is the distribution function (Bose-Einstein 
or Fermi-Dirac) and the integral is evaluated on the 
freeze-out surface E/ . In the SPH representation, Eq. 
|[5J| is written as 
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where the sum is over all SPH-particles on the freeze-out 
surface and n is the normal four-vector to this surface, 
given by oc (~dT/dr, -dT/dx, -dT/dy, -dT/drf). 



III. EQUATIONS OF STATE 

Let us consider here two sets of equations of state. 
The first one type-I (EOS-I) has a first-order phase tran- 
sition between an ideal gas of massless quarks (u, d and 
s) and gluons and an ideal hadron gas (baryon number 
is assumed to be zero), where the pressure, the energy 
density and the entropy density are given by: 



P q = 


—g q T 4 - B, 
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Here = 16 is an effective parameter (see [l7| for de- 
tails), g q — 95/2, and B is the bag model parameter. 

The second one, type-II (EOS-II), somewhat more re- 
alistic than the previous one, considers a first-order phase 
transition between a QGP and a hadronic resonance gas 
(baryon number is taking into account). In the QGP we 
consider an ideal gas of massless quarks (u, d, s) and 
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The hadronic phase is composed of resonances with mass 
below 2.5 GeV/c 2 , where volume correction is taken into 
account. The corrected- volume pressure, P ex , is written 
in function of the pressure for an ideal gas P ld : 



^^(T lMi )=^P/ d (T,/ii), 



(17) 



where fn = ^ — ViJ2jPj (T,jj,j), (V^: volume of the 
hadron i). Others thermodynamical quantities are given 
following the relation: 



l + ^VjnfiT,^)' 



(18) 



where \ represents n, e or s. 

In both equations of state, the transition temperature 
is assumed to be 160 McV. 



IV. INITIAL CONDITIONS AND RESULTS 

In order to study the behaviour of the effective tem- 
perature, we began with Landau- type initial conditions, 
where the fluid is at rest at t = fm/c, and the mat- 
ter is localized in a Lorentz contracted sphere. The ini- 
tial energy density is assumed to be constant and given 
by: £ = E/V, where E — C(v / s — 2m)A is the avail- 
able kinetical energy in the collision, and V — Jq 1 Vq — 
7 ( ^" 1 (4/3)7ri? 3 is the volume of the contracted incident nu- 
clei when superposed. Here £ is the inelasticity parame- 
ter and 70 = \fs/2m. These initial conditions (together 
with EOS-I) can reproduce both the particle multiplici- 
ties and the rapidity distributions fairly well (Ref. 17]). 
However it has an inconvenience that the initial energy 
density is too high. As a consequence, the time of ex- 
pansion becomes very large and transverse expansion is 
considerable. The result is a very large effective tem- 
perature and a non-exponential shape of the spectra in 
transverse momentum. 

A solution we found to this problem was to give the 
fluid an initial longitudinal-velocity distribution. We did 
not change the shape of the region where the fluid is 
formed, which remained a Lorentz contracted sphere. 
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FIG. 1: Energy dependence of the effective temperature for 
K + . The initial conditions are Landau- type ones, with initial 
longitudinal velocity. Triangles corresponde to Ty o =155 MeV. 
Experimental data are shown in square symbols. 



FIG. 2: Energy dependence of the effective temperature for 
K + . The initial conditions are NeXus-type ones. Triangles 
corresponde to T/ D =155 MeV. Experimental data are shown 
in square symbols. 



The initial longitudinal velocity is not constant, but pro- 
portional to z: 



u»{z) = ( 7 (z),0,0, 7 (z», 
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where z m = R/^o- As a consequence of this change, 
the initial energy density £o , which we took constant and 
obtained solving the equation 



E = 



+ p)j 2 (z) -p] dV, 



(22) 
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became much smaller and then the transverse expansion 
compatible with data, as shown in Table [i] and Fig. ^ 
In this simulation we consider freeze-out temperatures 
smoothly increasing with energy, until SPS domains. For 
RHIC energies, we show two cases: a decreasing freeze- 
out temperature and a increasing one (bold-faced types 
in Table [I] and triangle in Fig. QJ. 

However, since Sq became much smaller than the pre- 
vious case, so did the entropy density s an d the total 
multiplicity, although the rapidity distributions remained 
more or less the same in the shape. 

To reproduce all quantities (multiplicities, rapidity dis- 
tributions and transverse-momentum distributions) we 
have to take a large initial volume, with some appropri- 
ate longitudinal-velocity distribution. To do this, we use 
the NeXus event generator, which produces the event-by- 
event initial conditions at time r = 1 fm/c. Instead of 
using these fluctuating initial conditions, here we smooth 



TABLE I: Results obtained from Landau-type initial condi- 
tions with initial longitudinal velocity, vt is the transversal 
velocity averaged in the rapidity interval —0.5 < y < 0.5. eo 
is the initial energy density and To is the corresponding initial 
temperature. 





To 


eo 


Tfa 


VT 


T* 


(A-GeV) 


(MeV) 


(GeV/fm 3 ) 


(MeV) 




(MeV) 


2.7 


109 


0.097 


85 


0.240 


143 


3.3 


128 


0.185 


94 


0.255 


164 


3.8 


140 


0.263 


97 


0.280 


178 


4.3 


149 


0.341 


115 


0.155 


167 


4.9 


158 


0.431 


120 


0.153 


171 


8.8 


160 


1.038 


143 


0.146 


213 


12.3 


160 


1.322 


147 


0.147 


221 


17.3 


160 


1.536 


149 


0.139 


230 


130 


167 


1.896 


128 


0.292 


300 








155 


0.191 


273 


200 


168 


1.908 


125 


0.485 


392 








155 


0.272 


331 



them out by averaging 30 random events for each colli- 
sion energy. For this case we use the EOS-II. The results 
are shown in Table |n] and Fig. [21 Besides reproducing 
the effective temperature data quite well, it was shown 
that this initial conditions give good results also of mul- 
tiplicity and rapidity distributions |l8j . 
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TABLE II: Results obtained from NeXus-type initial condi- 
tions. The initial energy density eo and the initial tempera- 
ture To are given in the central point (x = y = rj = 0). vt 
is the transversal velocity averaged in the rapidity interval 
-0.5<y<0.5. 





To 


£0 


Tfo 


vt 


T* 


(A-GeV) 


(MeV) 


(GeV/fm 3 ) 


(MeV) 




(MeV) 


2.7 


98 


0.75 


85 


0.067 


92 


3.3 


128 


0.66 


94 


0.28 


155 


3.8 


131 


1.01 


97 


0.41 


192 


4.3 


135 


1.38 


115 


0.37 


212 


4.9 


140 


1.55 


120 


0.39 


225 


8.8 


198 


4.06 


143 


0.32 


234 


12.3 


248 


9.04 


147 


0.32 


231 


17.3 


265 


11.37 


149 


0.32 


228 


130 


279 


12.86 


128 


0.52 


280 








155 


0.34 


234 


200 


277 


12.48 


125 


0.56 


300 








155 


0.34 


231 



V. CONCLUSIONS 

We conclude that the initial conditions given in small 
volume like Landau-type discussed above are not appro- 
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priated to describe the multiplicity data, rapidity distri- 
butions and the transverse-momentum distributions at 
the same time. If we try to obtain the correct multi- 
plicities (together with the rapidity distributions), the 
effective temperature becomes too large (the first case 
discussed in Sec. IIVI where the fluid is at rest) . On the 
other side, if we try to get the correct effective temper- 
ature and the shape of rapidity distributions as we did, 
the multiplicities become small, unless an additional en- 
tropy is generated during the expansion (the second case 
discussed in Sec. IIV(I . It is clear that these conclusions 
about the initial conditions do not change if we replace 
EOS-I by EOS-II. 

Nexus type initial conditions, with spatially spread en- 
ergy and velocity distributions, reproduce well all the 
main characteristics of data, namely particle multiplici- 
ties, rapidity distributions and px spectra. In this case, 
the increase of T* in RHIC energy domain seems to be 
due mostly to the larger expansion in the hadronic phase 
than to that in the QGP phase. 
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